Hysteresis loops of magnetic thin films with perpendicular anisotropy 
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We model the magnetization of quasi two-dimensional systems with easy perpendicular (z-)axis 
anisotropy upon change of external magnetic field along z. The model is derived from the Landau- 
Lifshitz-Gilbert equation for magnetization evolution, written in closed form in terms of the z com- 
ponent of the magnetization only. The model includes-in addition to the external field-magnetic ex- 
change, dipolar interactions and structural disorder. The phase diagram in the disorder/interaction 
strength plane is presented, and the different qualitative regimes are analyzed. The results compare 
very well with observed experimental hysteresis loops and spatial magnetization patterns, as for 
instance for the case of Co-Pt multilayers. 
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I. INTRODUCTION 

Hysteresis is widespread in nature. Its origin and 
even what we mean by it may depend on the particular 
contexti Our concern here is that type of hysteresis ap- 
pearing due to the complex structure of the energy land- 
scape of the system studied, when some external control 
parameter is changed. In the paradigmatic description of 
the configuration space, the system is represented by a 
point evolving in a complicated energy landscape. If the 
available thermal energy is small compared to the heights 
of the typical barrier between local minima, the system 
can only evolve by moving down in energy and it will re- 
main in a local energy minimum. If some external control 
parameter is changed, the system adapts to the external 
conditions in order to reduce its energy. Upon changes 
of the external parameter it moves smoothly following 
the shift of the minimum in which it is located, until this 
minimum eventually becomes unstable. At this point the 
system rolls down in energy and finds another minimum 
in which it stays while this is stable. The evolution of 
the system in passing from one minimum that becomes 
unstable to another minimum implies energy dissipation. 
If the external field performs a cyclic variation, returning 
to its original value, the energy dissipated during the pro- 
cess can be calculated as the work done by the external 
field. If the external field B is conjugated to an internal 
variable M in the system, such that there is a term in the 
(differential) free energy of the system of the form MdB, 
then the evolution of the system in the AI -B plane is hys- 
teretic. Magnetic hysteresis is the prototypic example of 
this behavior, in which M is the system magnetization, 
and B is the external magnetic field. In this description 
of hysteresis the external field is assumed to change in- 
finitely slowly, in such a way that the system adapts in 
every moment to the actual value of the external field. 
However, this in general is not enough to guarantee a 
smooth evolution of the microscopic constituents of the 
system, and that is the origin of dissipation and at the 
end, of hysteresis. 

Although the qualitative origin of hysteresis is well un- 



derstood, a completely different issue is to reproduce in 
detail the observed hysteretic behavior of particular sys- 
tems. For magnetic materials this is an extremely im- 
portant issue, since hysteresis is at the base of magnetic 
information storage technology, and then its detailed un- 
derstanding has extreme practical importance. From an- 
other perspective, the different macroscopic character- 
istics of the hysteresis loops correspond to different mi- 
croscopic evolution of the spatial magnetization distribu- 
tion. The details of this correlation are worth to be stud- 
ied to gain insight into the fundamental physical princi- 
ples involved. 

In the present paper we model the hysteretic behav- 
ior of thin magnetic films with perpendicular anisotropy. 
In the next Section we discuss how to obtain a closed 
equation for the evolution of the perpendicular compo- 
nent of the magnetization, starting from the vectorial 
Landau-Lifshitz-Gilbert equation. The equation we ob- 
tain is qualitatively similar to some phenomenological 
proposals, but the derivation given here indicates clearly 
its possibilities and limitations. Section HI contains the 
results of the numerical simulation, and in Section IV 
we compare the results with experimental data in Co-Pt 
multilayers and conclude. 



II. DERIVATION OF THE MODEL 

Magnetism originates in quantum effects occurring at 
the atomic level. To model macroscopic hysteresis how- 
ever, it is impossible in practice to start from these ele- 
mental building blocks in order to arrive to the macro- 
scopic description. The micromagnetics approach}'^'' 
starts from an intermediate scale description, in which 
the magnetic moment of an elemental piece of material 
m is taken as the fundamental variable, m is a vector 
field, which is a function of the spatial coordinates and 
time, and it is assumed to satisfy at any position and 
time the constrain |mp ~ 1, where a rescaling of the 
amplitude has been assumed. A reasonably detailed de- 
scription of magnetization evolution is provided by the 
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phenomenological Landau-Lifshitz-Gilbert equatioiifi*2i^ 

that can be written as 



9m 



= am X B — 6m x (m x B) 



(1) 



where B is the local effective magnetic field. The first 
term on the right of Eq. describes a precessional 
evolution of m around B, and the second term is a phe- 
nomenological damping term that tends to align m with 
B. Note that from this equation we obtain immedi- 
ately that d\m\'^/dt = 0, i.e, it automatically maintains 
the norm of m as fixed. In equilibrium {dm/dt = 0), 
this equation reduces to Brown's equation,^'^ namely 
m X B = 0, stating that at any point the local mag- 
netization has to be aligned with the local field. The 
local field B includes contributions from the externally 
applied field, from the interaction with all the rest of the 
magnetic moments in the sample, and from structural 
anisotropics in the sample. 

Equation is a sufficiently detailed starting point to 
obtain the time and space evolution of the magnetiza- 
tion. Unfortunately, it is frequently a too much detailed 
starting point. In fact, except in those simple cases in 
which an analytical approach is possible, the numerical 
solution of Eq. ^ is difficult because the orientation of 
the local field changes in general with time. 

For the problem we want to study we can obtain from 
Eq. a numerically tractable simplified formulation. 
The problem is that of the magnetization evolution of 
(quasi) two-dimensional systems with perpendicular (z) 
easy axis anisotropy, in the presence of an external mag- 
netic field along z. This problem has attracted a lot of 
attention, in particular because it constitutes a promis- 
ing configuration for high density magnetic recording."* 
The perpendicular easy axis anisotropy induces that in 
most of the sample the magnetization points in the plus 
or minus z direction, defining magnetization domains. 
Only within the 'domain walls' separating magnetic do- 
mains with opposite magnetization, the local magneti- 
zation cannot be considered to point along z. In the 
assumption that the extent of the system occupied by do- 
main walls is small compared to the domains themselves, 
we can derive an effective model from Eq. ^ taking into 
account only the z component of the magnetization. 

To derive this model, we first of all need to specify in 
more detail what the contributions to the local field B 
are. We will consider the following four contributions to 
B: 

f m(r') 

B = hoi + Am,z - -f dr' '—L^+Bexch (2) 

The first term is an external field applied along z. Second 
term is the lowest order contribution from the assumed 
easy axis anisotropy of the sample. Third term is the 
dipolar (or stray) field generated by all the rest of the 
sample. The last term accounts for the exchange field. 
It typically consists of a sum of the magnetization values 



on neighbor sites, and then it introduces an explicit de- 
pendence on the symmetry of the underlying lattice. In 
the simplest case in which we want to get rid of lattice 
effects to concentrate on the intrinsic dynamics of the 
model, it is convenient to have a symmetric form for the 
exchange term. By noticing that if we add to B a term 
proportional to the local magnetization it does not affect 
the equation of motion (since only B x m appears), the 
exchange term can be written in the form 



B 



exch 



(3) 



where J defines the exchange constant. Note that the 
fact that we take |m| = 1 is incorporated in the values of 
A, 7 and J. 

To obtain a closed equation for the evolution of niz, 
we notice that under the assumptions made, the dipolar 
field points to a good approximation along the z direc- 
tion everywhere (since non-z components are originated 
in the field generated by the domain walls, which is much 
smaller). Then the first three terms of Eq. Q point in 
the z direction. We separate explicitly this part from the 
exchange part: 



B ~ {ho + Am,,-j I rfr'^^^^ ] z + JV^m 



hz + Bexch 



(4) 



from these two terms we get two contributions to the z 
component of Eq. 



dr 



dt 



= b{l—ml)h— am x B 



Although within the domains B^xch 
rection, it certainly does not within domain walls 



exch\:,-b m X (m X Bexch)]:, 

(5) 

points in the z di- 
The 



B, 



term will infiuence the detailed structure of the 



domain walls. Since we are not interested in the in- 
ternal structure of the walls, we can proceed by replac- 
ing Bexch by some phenomenological expression depen- 
dent on niz only. One way of obtaining such an expres- 
sion is the following. Consider a domain wall extending 
along the y axis. If we imagine that the internal struc- 
ture of the wall is that known as a Bloch's wall, then 
TTix = 0. On the contrary, if we think that the wall is 
a Neel's wall we have niy = Oj^, In both cases magneti- 
zation can be described in terms of the angle 9 between 
local magnetization and the z axis, in such a way that 
niz = cos{9), and my = ±sin(6'), nix = for a Bloch's 
wall, or nix = ±sin(6'), my = for a Neel's wall. In 
any of these cases the second term in the r.h.s. of Eq. 
© still vanishes, and the last term takes the simple form 
6Jsin(6')^, that can be re-expressed in terms of and 
reinserted into (O to give 



dm;, 
dt 



= b 



(1 - ml)h + J 



dx^ 



J 



( dm^ 



1 — ml \ dx 



(6) 
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Once in this form, the equation can be generahzed to 
describe a system with domain walls oriented in any di- 
rection, in the form 



dr 



dt 



= b 



(1 - mj) ha + Am.;, - 7 / dr' 



+J\7^m;, + J 



1 



r(VTOz.VTO^) 



(7) 



where we explicitly introduced the form of h from Eq. 

As we explained above, this equation will be mean- 
ingful only in cases in which the magnetization distribu- 
tion consists of domains of positive an negative magne- 
tization, separated by domain walls, and such that the 
thickness of the domain walls is small compared to the 
typical width of the domains. We first present some sim- 
ple examples of the use of this equation, and then the 
results of numerical simulations to study the form of hys- 
teresis loops of two dimensional garnets. 

The first simple example is that of a single, isolated 
spin in the presence of an external field hp which is fixed 
in orientation (along z), although its amplitude may de- 
pend on time ho = ho{t)z. The only term in the r.h.s. of 
Eq. (O that survives in this case is the first one, and the 
evolution equation is: 



dm^ 



= h{l-ml)h„{t) 



(8) 



that can be immediately integrated to give 



'niz{t) = tanh I b I ho{t)dt 



(9) 



The particular case in which ho is constant describes the 
spin flip process of an isolated spin in a constant field. 

The second example is the stationary structure of 
a straight domain wall in the presence of structural 
anisotropy, with no dipolar or external field. The equa- 
tion to be solved becomes 



A{1 — m1)mz = J 



d mz 



1 



dm^ 
dx 



(10) 



which in terms of the variable 9 (such that cos(6') — m^) 
is written as 



= — sm[&) cos(t') 

and from here the well known solution 

A 
1 



mz = cos{9) = tanh x 



(11) 



(12) 



is obtained. We see from here how the thickness of the 
wall is proportional to the ratio J/ A between the ex- 
change and anisotropy constants. It can be verified that 
the same spatial dependence but with a domain wall 



thickness which is doubled, is a solution of Eq. (|10|l if we 
drop the last term, quadratic in dmz/dx. Below we will 
return to this point. 

Before moving to the simulations of two dimensional 
systems, we will compare the present model with oth- 
ers that have been proposed to describe similar prob- 
lems. Starting from the basic Ising model of discrete ±1 
spins, Sethna et al.^ have introduced the random field 
Ising model (RFIM). This is an important model system, 
where the effect of disorder on the existence or not of an 
abrupt nucleation step was analyzed in detail. Other 
people have considered the effect of adding dipolar in- 
teraction to the RFIM. The dipolar RFIMSi contains in 
principle all physical ingredients necessary to provide a 
good description of uniaxial anisotropy systems (not nec- 
essary thin films) in the presence of magnetic field along 
the easy direction. However, it suffers from a technical 
problem, namely that since the fundamental variable is 
discrete, domain walls are forced to be one lattice pa- 
rameter thick. This produces domain walls that are ar- 
tificially pinned to the numerical lattice, and this makes 
realistic simulations very difficult. An alternative to over- 
come this problem consist in using a continuum variable 
0, which instead of taking only the values ±1, can have 
intermediate values, at the cost of paying some energy. 
We obtain in this way what can be called the "dipolar 
phi-fourth" model, described by the equatior '''^'^ 



d<j)[r) 
dt 



^ho-A{-cl) + c^^)--/ / dr' 



r'|3) 



(13) 

The scalar variable (p (which may be assimilated to the 
z component of the local magnetization) evolves in a lo- 
cal two-well potential which mimics the existence of two 
privileged values (±1) for (j). External field, dipolar and 
exchange contributions are included. We see a general 
similarity between this equation and our Equation l|7|l. 
However, a few differences need to be mentioned. The 
most obvious one is that in Eq. (|13|l (p is not restricted 
to satisfy |(/)| < 1 and in fact, larger values are obtained, 
in particular if the external field is strong enough. In 
Eq. {TJ, the restriction jm^l < 1 is guaranteed by the 
(1 — ml) factor. A second difference is in the structure of 
the gradient terms: the last, non-linear term in Eq. Q 
does not appear in Eq. H13|l . Although in some cases this 
term may have some physical importance, we have seen 
that for a simple, static domain wall, it only renormalizes 
its thickness. In view of the difficulties that the presence 
of this term would add to the simulation (making them 
much more time consuming because of its non-linear na- 
ture) we have decided to drop it in the simulations that 
are presented below, and assume that its main effect has 
been incorporated in an appropriate renormalization of 
the amplitude of the V^m^ term. This simplifying as- 
sumption is justified also since we do not intend to model 
details of the interior of the domain walls, and then all 
that we need is a term that produces some energy cost for 
the domain wall. The V^m^ term suffices in this respect. 
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III. RESULTS 

The model simulated is thus given by the equation 
dm 



2.1 



J/Ao 1 



dt 



r'|3 



(14) 

which is Eq. {T)) in the absence of the last term, where 
time has been rescaled to absorb the constant b, and 
where a spatial dependence of the coefficient A has been 
made explicit. This dependence will allow us to simulate 
structural disorder in the sample. A single, isolated spin, 
is in equilibrium if it points in the direction of the local 
field. We see that rriz — +1 {niz = —1) is an equilibrium 
configuration of the spin if /io+ A > {ho — A < 0). Then 
A represents the value of the local coercive field that is 
necessary to apply in order to invert the orientation of 
the spin. The function ^(r) will be chosen in the form 



Air) ^ Aoil + Dr^r)) 



(15) 



where 77 is a spatially random and uncorrelated Gaussian 
variable, with zero mean and unitary variance {rj is cut 
off at large negative values, namely rj > — 1/D, to guar- 
antee A > 0). D controls the overall intensity of disorder. 
The parameters of the model are thus Aq, D, 7 and J, in 
addition to the external field hg. As explained in [9], the 
ratio between J and 7 can be adjusted by appropriately 
rescaling the spatial scale in the model. In the simula- 
tions below we choose 7 = 0.095 J, and give the results 
in term of ho/Ao, J/Aq, and D. The mesh parameter is 
taken as the unit of length. 

The results presented below correspond to zero tem- 
perature. Thus the evolution of the system is driven 
exclusively by the tendency to minimize the energy, as 
described in the introductory part. We will see that even 
this T = case is very rich, and provides different macro- 
scopic forms of the hysteresis curve when the interaction 
strength between spins and the amount of structural dis- 
order are changed. The results obtained compare very 
well with experimental results. 

The main result we are going to present is the disorder- 
exchange/anisotropy ratio (D-J/Aq) phase diagram of 
the model, and its description. The phase diagram is 
presented in Fig. ^ The form of the hysteresis loop at 
different positions is shown in Fig. |21 

There are three main qualitative regions in the phase 
diagram. For large J/Aq, (2.1 < J/Aq) the system shows 
no hysteresis at all. A closer examination of the magne- 
tization distribution shows that in this region the inter- 
action is so strong that effectively, the two minima of the 
local magnetization disappear, and tUz remains almost 
uniform in the sample, changing smoothly as a function 
of the applied field. This regime is certainly outside the 
region in which the model is justified, and then we will 
not consider it. 

More interesting to us is the rest of the phase diagram, 
where hysteresis occur, and where (from the examination 
of magnetization distribution) we observe well defined 
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FIG. 1: Disorder-exchange/anisotropy ratio (D — J/Aq) phase 
diagram for the model described by Eq. 11411 . in a system of 
size 384 x 384, with 7 = 0.095 J. Numbers indicate the 
points for which the hysteresis loops are shown in the next 
figure. The continuous line is the main feature of the diagram, 
separating regions in which an abrupt nucleation event occurs 
(to the left of this line) or not (to the right). In addition, the 
nucleation event produces complete magnetization reversal to 
the left of the dotted line, but only partial magnetization 
reversal between this curve and the continuous line. 



domains with magnetization ±1. That is the part plot- 
ted in Fig. 121 A mayor distinction can be made within 
this region. It is related to the existence or not of a finite 
magnetization jump during field cycling. The (numeri- 
cally determined) regions in which a finite jump exists 
or not are separated by the continuous line in Fig. ^ 
In addition, if a jump exists, it can lead directly to the 
fully inverted magnetization state, or to a partial inver- 
sion, that requires a further change of the external field 
to completely revert the magnetization. This two pos- 
sibilities are separated by the dotted line in Fig. ^ A 
magnetization jump indicates an instability in the sys- 
tem, in which a finite fraction of the spins change their 
orientations, upon an infinitesimal change of the exter- 
nal field. This avalanche process in magnetic systems 
was originally analyzed in the random field Ising model 
(RFIM).^ It was shown that if disorder in the system is 
lower than some critical amount, then there is a magne- 
tization jump, whereas there is no jump for higher disor- 
der. The same result is qualitatively obtained here, the 
amplitude of the magnetization jump vanishes continu- 
ously when we approach the continuous line from the left 
in Fig. n This line is a critical line, in the same sense 
used in the RFIM. We expect that the vanishing of the 
abrupt jump follows a power law close to the critical line. 
However, we have not attempted to determine critical ex- 
ponents since in the presence of the dipolar interaction 
we cannot go to system sizes large enough to get good 
statistics. An interesting question that remains open is 
whether the critical exponents are the same all along the 
critical line or not. 

The spatial distribution of the magnetization shows 
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FIG. 2: The hysteresis loops at different locations of the 
phase diagram shown in the preceding figure. 

In the vertical axis we plot the magnetization m^, which 
changes between ±1 as the field is changed. 

characteristic features for different forms of the hystere- 
sis loops. Examples are shown in Figs. |3 01 |5l |H| and 
[7| corresponding to hysteresis loops (1), (4), (9), (7), 
and (10) in Fig El Fig. corresponds to a case where 
there is an abrupt nucleation, leading to a partially in- 
verted magnetization state. Panel (a) in Fig. (31 is a non- 
equilibrium configuration after the nucleation. The final 
state at this field corresponds to a configuration of mag- 
netization ~ 0.3. Further decrease of the field (panels (c) 
and (d)) is necessary to completely invert the magneti- 
zation. In the presence of strong disorder, for approxi- 
mately the same value of J/Aq (Fig. 2} the nucleation 
step is smoothed, and the whole evolution is continuous. 




(c) m = -0.5 (d) m = -0.835 

FIG. 3: Four snapshots of the spatial magnetization distribu- 
tion corresponding to hysteresis loop (1) in Fig|51 The values 
of the magnetization at which the snapshot were taken are 
indicated. Note that the first panel is an unstable configu- 
ration, which is spontaneously evolving towards a state with 
m ~ 0.3. 



Here disorder plays a fundamental role in preventing the 
spontaneous growth of the nucleated domains, which in 
this case remain pinned by the structural disorder, and 
need a further decrease of the field to be able to grow. 

Fig. [31 displays a case in which the nucleation events 
takes the sample completely in the oppositely magnetized 
sated. Note how in this case, the domain that nucleates 
is a rather featureless, more or less circular objects, and 
grows to take over the whole sample. This case to be 
compared with Fig. (31 in which the domains nucleated 
maintain a striped internal structure. An intermediate 
case is that of Fig. El Here the inverted magnetization 
bubble almost invades the whole sample at once. How- 
ever, a small fraction of domain with the original ori- 
entation remains. This is reflected in the form of the 
magnetization loop, which in fact shows an abrupt jump 
of the magnetization from -1-1 to about -0.5. 

The last case (Fig. [TJ corresponds to a much weakly 
interacting sample (lower value of J/Aq). However, the 
disorder is sufficiently large to avoid a sudden nucleation 
step. We see a very different configurations of the do- 
mains in this case, which do not show the typical striped 
pattern of the other cases (originated in the dipolar inter- 
action) but patterns mainly dominated by the disorder. 
This case, in which the dipolar interaction play a minor 
role, is qualitatively similar to that studied in the RFIM. 

One striking characteristic of the phase diagram we are 
presenting is the reentrance of the sector with continuous 
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(c)m = -0.33 (d)m = -0.67 (c)m = -0.5 (d) m = -0.835 



FIG. 4: Same as Fig. 13 for hysteresis loop (4) in Fig El FIG. 6: Same as Fig. |21 for hysteresis loop (7) in Fig |5] 

First two panels show the spontaneously growing nucleated 
domain. 




(c) m = -0.5 (d) m = -0.835 

FIG. 5: Same as Fig. for hysteresis loop (9) in Fig |5| 
None of these configurations is stable, they are spontaneously 

evolving to the completely inverted magnetization state. FIG. 7: Same as Fig. El for hysteresis loop (10) in FigEl 



hysteresis. This reentrance tells that if D is not too high, 
samples displaying continuous hysteresis can be classified 
in two well different groups: weakly interacting ones for 
low J/Aq, and strongly interacting ones for large J/Aq. 
For strongly interacting samples, the width of the hys- 
teresis loop (measured as the distance between the two 



branches, at a fixed magnetization) can be much smaller 
than the change of magnetic field necessary for complete 
magnetization inversion (see for instance loop (3) in Fig. 

In these samples the coercive field is much lower 
than the saturation field, and remmanent magnetization 
is much lower than saturation magnetization. 
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In weakly interacting samples the range of external 
field on which the complete magnetization inversion is 
produced is small compared with the width of the hys- 
teresis loop (as for instance in loop (10) in Fig. |2J). As a 
consequence, these samples posses a coercive field which 
is of the order of the saturation field. In addition, the 
remmanent magnetization is essentially the saturation 
magnetization. 

The case of what we have called 'weakly interacting 
samples' allows yet a further distinction, which however 
cannot be inferred from the results shown up to now 
only. In fact, weakly interacting samples are examples 
of 'permanent magnets'. A usual distinction among per- 
manent magnets is that between nucleation type magnets 
an pinning type magnets.^ They are phenomenologically 
distinguished by the form of the virgin magnetization 
curve. Pinning type magnets display a small initial per- 
meability at the demagnetized state, whereas for nucle- 
ation type magnets the initial permeability is very large. 
The two different behaviors can be obtained with the 
present model. An example is shown in Fig. |S1 (for 
simplicity in the simulation, we obtain a demagnetized 
system by quenching from a hight temperature configu- 
ration, instead of the annealing process usually invoked). 
Pinning type magnets are modeled through a very low 
value of the interaction. The system behaves essentially 
as a collection of isolated magnetic moments, that switch 
in the presence of the external field at the particular value 
of the local coercive field. The transition is seen to be 
smoothed because of the disorder in the coercive field. 
For larger values of the interaction strength (but still 
in the region corresponding to 'weakly interacting sam- 
ples') we can see how the initial permeability increases 
noticeable. Here, domain walls (which are abundant in 
the demagnetized state) have some freedom to move in 
the presence of the external field, and this produces a 
high permeability. However, once fully magnetized and 
upon inversion of the external field, domains with op- 
posite magnetization have to nucleate, and this requires 
much larger fields. Note how the difference in the two vir- 
gin curves occur in two samples that have an extremely 
similar appearance of the full hysteresis loop. 



IV. COMPARISON WITH EXPERIMENTS AND 
CONCLUSIONS 

There has been in recent years much interest in the 
behavior of thin multilayer magnetic films^S*iiiSi mainly 
as potential high density magnetic storage media. Par- 
ticular attention is being paid to Co-Pt multilayers. The 
influence of disorder on these material has only recently 
started to be studied systematically. An example of the 
evolution of hysteresis loops upon change of the struc- 
tural disordered in the sample is contained in [12]. There, 
Co-Pt multilayers were grown under different values of 
argon sputtering pressure, what allows for the introduc- 
tion of disorder in a controlled way. We can compare the 




-2.0 -1.0 1.0 2.0 

FIG. 8: Hysteresis loop and virgin curve, obtained from a 
demagnetized sate, showing two different qualitative cases: a 
'pinning dominated' magnet (a) in which the initial perme- 
ability of the sample is small, and a 'nucleation dominated' 
magnet (b), with a large initial permeability. Note that both 
cases correspond to what has been called "weakly interacting 
samples" . 



plots in Fig. 1 of Rcf. [12] with the hysteresis loops for 
the present model for a value J/Aq ~ 2 (see panels (1) 
to (5) in Fig. |2Jl in which the disordered is progressively 
increased. The qualitative agreement between the model 
and the experiment is astonishingly good, taking into ac- 
count the simplifications in the model. Comparison of 
our results with images from magnetic scattering tech- 
niques (Fig. 2 in Ref. [12], for instance) reveals also that 
even the real space configurations are quite realistic. We 
see that the samples analyzed in [11] and [12] fit very well 
in the strong interacting region of our phase diagram. 

The relevance of the present results to experiments is 
twofold. On one side, it provides insight on what are the 
main physical ingredients that have to be considered to 
obtain a thorough description of the phenomenon. On 
the other side, it it also of importance in the following 
sense. In many cases, experimental information of the 
real space patterns is inferred from observation of the 
X-ray diffraction patterns, in which only the amplitude 
information is conserved. Then the conclusions drawn 
from the experiments depend on the possibility of ex- 
tracting information of the real space patterns from the 
X-ray patterns, and this is not a trivial issue. The present 
model, providing directly the configuration in real space 
(which can of course be transformed to get the X-ray 
patterns) is an ideal benchmark in which the X-ray re- 
construction techniques can be tested. A more detailed 
statistical comparison between simulated and measured 
patterns along these lines will be published elsewhereii^ 
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